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A numerical renormalization group (NRG) investigation of the one-centre t — t' Kondo problem 
is performed for the square lattice with account of logarithmic Van Hove singularities (VHS) in the 
electron density of states. The magnetic susceptibility, entropy and specific heat are calculated. The 
temperature dependences of the thermodynamic properties in the presence of VHS turn out to be 
non-trivial. For finite t' inverse logarithm of the corresponding Kondo temperature Tk demonstrates 
a crossover from the square-root to standard linear dependence on the s — d exchange coupling. The 
low-temperature behavior of magnetic susceptibility and linear specific heat are investigated, and 
the Wilson ratio is obtained. For t' Q the Fermi-liquid behavior is broken. 
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I. INTRODUCTION 

The Kondo effect is extensively studied starting from pioneering works by Kondo [l] who explained the minimum 
of resistance in metallic alloys owing to resonance s — d scattering of conduction electrons by magnetic impurities. 
The solution of the Kondo problem ^2n^] is a very beautiful chapter in the history of modern theoretical physics. 

The Kondo phenomenon is a key for explaining the behavior of heavy- fermion compounds and Kondo lattices Q , 
non- Fermi- liquid (NFL) systems [3, @], anomalous electronic properties of metallic glasses [9], quantum dots [l0|, and 
other systems. The Kondo anomalies are also studied in systems of reduced dimensionality where they have a number 
of experimental peculiarities [ll| . It is evident that the Kondo effect in such systems has a number of peculiar features 
from the theoretical point of view too. 

The Kondo effect owing to Cu^"*" spins in the Cu02 planes is observed in layered n-type cuprates (La,Ce)2Cu04, 
(Pr,Ce)2Cu4, and (Nd,Ce)2Cu04-5 

Last time, the Kondo effect in graphene (truly two-dimensional system with a peculiar electron spectrum) is 
discussed [l3l - [l5t . Because of the pseudogap in the spectrum, the Kondo effect for the undoped graphene exists under 
restricted conditions only, but for a doped substance the Kondo phase is present for all parameters [I^l ■ 

We can mention also some layered /-systems where experimental investigations and first-principle band calculations 
demonstrate existence of two-dimensional features in electron properties. Here belong the compounds CeCoIus (where 
the layers Cein are present) [lB|, CeCoGe2 



, _ , , - . 17|, CePt2ln7 [18|, CeRhlng, Ce2RhIn8 [l% UCoo.5Sb2 (where two- 

dimensional weak localization is observed) ^ | . 

The layered Kondo lattice model was proposed for quan tum critical beta-YbAlB4 where two-dimensional boron 
layers are Kondo coupled via interlayer Yb moments [21[. CeRuPO seems to be one of the rare examples of a 
ferromagnetic Kondo lattice where LSDA+U calculations evidence a quasi-two-dimensional electronic band structure, 
reflecting a strong covalcnt bonding within the CeO and RuP layers and a weak ioniclike bonding between the layers 

M- 

The above /-systems demonstrate often both local-moment and itinerant-electron features. Large linear specific 
heat and NFL behavior is observed also in some d-systems including layered ruthenates Sr2Ru04 [23| and Sr3Ru207 
[23 |. Besides correlation effects, anomalies of electron properties in such systems are owing yo the presence of Van 
Hove singularities near the Fermi level. 

In the present work we treat the one-centre Kondo problem with the singular electron density of states. Earlier 
this problem was considered by Gogolin [25] who used a "fast parquet" approach. In fact, such complicated methods 
are somewhat ambiguous, and the numerical renormalization group (NRG) results of paper [2^ do not agree with 
the results of Ref.[2^. Therefore we start in Sect. 2 from the standard perturbation theory and also apply the "poor 
man scaling" approach by Anderson [27| . 

In Sect. 3 we apply to the problem the NRG method, the technical details for our case being considered in Appendix. 
The simple perturbation results for the Kondo temperature agree with NRG much better than the parquet results 
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[25| . The physical quantities in the presence of the logarithmic singularity near the Fermi level are investigated. 
We calculate the magnetic susceptibility and specific heat, in particular at low temperatures, discuss the problem if 
universal behavior and calculate the Wilson ratio. 



II. THE KONDO MODEL WITH VAN HOVE SINGULARITIES 



We use the Hamiltonian of the one-centre s — d{f) exchange (Kondo) model 



(1) 



kCT 



where £k is the band energy, S are spin operators, a are the Pauli matrices, in the case of contact coupling Jkk' — J/Ns 
where J is the s — d{f) exchange parameter, Ns is the number of lattice sites. 

Kondo f]\ found that high-order perturbation contributions to various physical properties contain logarithmically 
divergent corrections. As demonstrated further investigations of the Kondo problem, there occurs a pole at the 
boundary of strong coupling region, which is called the Kondo temperature (in fact, this is a crossover scale). For a 
smooth density of states p{E) this quantity is estimated as 



Tk cc D exp ; 



1 



2 Jp(0) ' 

where D is the half-bandwidth. We treat the case of logarithmically divergent bare density of electron states 

D 



(2) 



p{E) = Aln 



B\E + A\ 



(the energy is referred to the Fermi level, the constants A, B and A are determined by the band spectrum). The 
logarithmic divergence in p{E) is typical for the two-dimensional case (in particular, for the layered ruthenates). 
However, similar strong Van Hove singularities can occur also in some three-dimensional systems like Pd alloys and 
weak itinerant ferromagnets ZrZn2 and TiBr2 [28, 29]. 

First we consider perturbation expansion for the resistivity, following to the original approach by Kondo [1,]. We 
write down the inverse transport relaxation time with the Kondo correction 



-\E)=T^HE)[l+AJg{EM- 



Here 



g{E,T)^Y. 



1/2 ~ 
E-ei, 



= dE'p{E') 



,^l/2-/(i?') 



E-E' 



(3) 
(4) 



with rik = /(^k) is the Fermi function. After integration by part, in the case A = we obtain to logarithmic accuracy 
for the resistivity 



1 - 2AJ dE' i - 



df{E'] 
dE' 



In^ 



D 
E' 



so that 



9 D 

T 



R 



(0) 



J^S{S+l)A\n^ 



Applying the Abrikosov-Suhl summation (see Q) we get 



Iisd = RfJ (1 + J^ln 



D 
T 

which yields a non-standard expression for the Kondo temperature, 

Tx — D exp 





1 


1/2- 




AJ 





(5) 



(6) 



(7) 



(8) 
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We calculate also the Kondo corrections to the static impurity magnetic susceptibility by generalizing consideration 
of Ref.fll to the case of singular density of states. Expanding to second order in J we derive (cf. [ll, fflsoj) 



xiT) 



S{S + l) 
ST 



kk' 



^T-k(l - ?^k' 
(Ck -e-k' f 



where 



D 



max(A - T, T) 



(9) 



(10) 



is the Pauli susceptibility of non-interacting conduction electrons with the singular density of states, which is shown 
in FigH] For A ^ this quantity has a maximum at T ~ A/2. Such a behavior is typical for the case where a 
density-of-states peak is present near the Fermi level jsij . 
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FIG. 1: The temperature dependence of non-interacting magnetic susceptibility for conduction electrons X^'^^T) at A = 0, 
0.02, 0.04, 0.08 



Performing integration and summation of the series of logarithmic terms we have 



X{T) 



S[S + 1) 
3T 



2Jx(0) 



1 + Jx^°'> ln{D/T) 



(11) 



which yields at A = the same result for the Kondo temperature ([8|) 

The correction to magnetic impurity entropy can be written down in an analogous way to obtain [l|, 

^ 5(5 + 1) [2Jx(°)l' 



(12) 



Then we have for impurity specific heat 



amp(T) = T 



dT 



7r25(5+ 1) 



[1 + Jx("' HD/T)] 



2JX 



(0) 



dT 



.(0) 



(13) 



The first term in the brackets yields the same structure as in the case of a smooth density of states (where x^^^ — P ^-^d 
the singular contribution occurs in the fifth order in J only) , and the second term is owing to logarithmic singularity 
in x^*^^. Being of the third order in J, the latter term can dominate. Beside that, it can change its sign and become 
negative; as we shall see below in Sect. 3, this is important for the Cimp(r) behavior. 
Of course, the above expressions are applicable for T > Tk only. 

To perform a more formal consideration, we can apply the "poor man scaling" approach [13]. This treats the 
dependence of effective (renormalized) model parameters on the cutoff parameter C which occurs at picking out 
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the singular contributions from the Kondo corrections to the effective coupUng Jef{C) (with Je/(— D) = J). The 
second-order singular correction to Jef{C) can be obtained in the form [s^ 



(14) 



Picking out in the sum the contribution of intermediate electron states near the Fermi level with C < < C + SC 
we obtain 



SJefiC) = 2Aj2ln- 

Replacing in the right-hand part J — > Jef (C) we obtain 

d 1 2A 



D 5C 
C + A\~C 



dCJefiC) 



In- 



D 



C \C + A\ 



Solving this equation for Tk ^ A we get 



' ^1 + Aln^ 



Then we obtain from the condition l/Je/(Tx) = again the result ([5]). 

One can see that the expression (|8]) is different from the corresponding parquet result [2£ 



Tk — D exp 
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(15) 



(16) 



(17) 



(18) 



However, the NRG calculations [23 confirms the perturbation expression ([SJ rather than ([T5| (see the discussion 
below). The corresponding problems of the parquet approximation in the Hubbard model are discussed in the works 

fH. 



III. RESULTS OF NUMERICAL CALCULATIONS 



We consider 5=1/2 Kondo model for the square lattice with the spectrum 

£k — 2t(cos kx + cos ky) + 4i'(cos k^ cos ky -^- 1) 
The corresponding density of states is 

1 



P{E) 



27rVi2 + i^t' - 4t'2 



K 



1 



In 



t^ + Et'+4t'^ 27rVi2 - 4i'2 \E + 8t' 



(19) 



where K{E) is the complete elliptic integral of the first kind. In the numerical calculations, the Fermi level is supposed 
to be located in the band centre, so that the band is determined by \E\ < D = A\t\. The distance from the Van Hove 
singularity to the Fermi level is A = 8t' . When presenting numerical results, we put D = 1. 
For t' = we derive 



piE) 



2 , AD 
■ In ■ 



TT^D E 



(20) 



so that, according to 



Tk — D exp 







1/2- 




2J 





(21) 



The details of NRG calculations are discussed in Appendix. 
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FIG. 2; The dependence Tk{J) for t' = 0. Circles are NRG results, solid line corresponds to the Nagaoka-Suhl formula (|22p . 
dashed line to Eq.((2l| and dotted one to Eq. (|18|) (without any fitted preexponential constants) 



The Kondo temperature in NRG calculations is determined from the temperature dependence of impurity magnetic 
susceptibility Ximp 

(T) by using condition rKXimp(Tx) = 0.0701 that is standard in the NRG method 3. 

Figs mo show the dependence of the Kondo temperature on the bare coupling | J| in the logarithmic scale. One can 
see that for finite t' there occurs a crossover with decreasing | J| from the square-root dependence ([5]) to the standard 
Kondo behavior ([2]). This crossover is qualitatively described by the lowest-order scaling. 

It should be noted that the Kondo phenomenon should be invariant under particle-hole transformation (ek — ^ ^Ck, 
i.e., p{e) — p(— e)). However, in the case of a non-symmetric density of states using the formula 1 = 2Jg(TK,0) 
(where g{E,T) is determined by yields slightly different results for the peak above and below the Fermi level. 
Therefore we use the Nagaoka-Suhl formula (see [l|) 

l = 2Jg{0,TK) (22) 

which works somewhat better at large J, but slightly worse at small J. (In fact, all such approximations yield nearly 
the same result at small J.) 




FIG. 3: Fitting of the dependences Tk{J) for t' — 0, 0.005, 0.01, 0.03. Circles, top-up triangles, squares and top-down triangles 
are corresponding NRG results, solid lines (from above to below) correspond to the Nagaoka-Suhl formula (I22|l . dashed lines 
(from above to below) show the two-loop result (I23p with fitted constants k 



According to p^ . for A <^ Tk we have InT^ ^ l/\/PT ^-iid for A ^ Tk we come to the standard linear behavior 
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IhTr- ~ with enhanced p(0). For intermediate J, a crossover takes place. 

The lowest-order scaling describes satisfactorily the numerical data at large | J|. However, it is insufficient to fit the 
numerical results at small |J|. Therefore we use the two- loop scaling result 



Tk = fcD|2Jp|i/2gxp[i/2Jp] 



(23) 



This expression implies that only small vicinity of the Fermi surface with nearly constant p{E) ~ p works. This 
assumption becomes not valid with increasing \ J\ since the whole logarithmic peak becomes important. It should be 
noted that, as follows from the structure of scaling equations (see Ref.Q), the factor of \2J p\^/'^ does not occur in the 
expression for Tk at A = 0. 

The fitted constant k should be determined by the whole form of the function p{E). At small t' , k is small since 
the quantity p{0) is large owing to peaks. On the contrary, for larger t' — 0.03 — 0.05, p(0) becomes small. 

Since Tk is high in our case of singular density of states, consideration of the situation with A ~ Tk is quite 
realistic. We can see that appreciable deviations owing to the singularity occur even if A is not too small and exceeds 
Tk- 

The magnetic susceptibility owing to impurity can be expressed as a difference of magnetic susceptibilities of the 
whole system and the system without impurity: 



X{T) = Ximp(r) = Xtot(r) - Xband(r), 



(24) 



where xtot is the total magnetic susceptibility, and Xband = 2Afsx'''^' is the susceptibility of non-interacting band 
electrons (for two spin projections). Apart from the susceptibility ([Ml), the so-called local magnetic susceptibility xioc 
is frequently introduced as well: 



l/T 



Xloc(T) 



(5,(r)5,)dr 



(25) 



This is the susceptibility of a single impurity in a magnetic field that acts locally, only on this impurity; its magnitude, 
therefore, can hardly be measured experimentally, in contrast to Ximp- 

In principle, Ximp E^nd xioc can behave quite differently. This possibility was mentioned in Ref.fsJ], where the 
reason for the difference is related to the energy dependence p{E) and is asserted that this disappears for a flat band 
of half-width D in the limit of 13 — > oo. 

The impurity entropy and specific heat are defined in a similar way. 

Figs mil] show the temperature dependence of magnetic susceptibility for different t' values. At high temperatures, 
x{T) obeys the Curie law, and at low temperatures it demonstrates the Pauli behavior (except for the case A — )• 0). 



E 




FIG. 4: The temperature dependence rximp(r) for t' = and J — —0.2, —0.1, —0.07, —0.05 (lines from below to above) 



For the standard flat-band case we have Tx{T) — 4>{T/Tk), so that the curves Tx{T) are universal: a change in 
J results in a change of Tk only. In our situation, such a simple universality does not hold. In particular, for t' = 
this fact (illustrated by FigH]) was demonstrated in Ref . [26l| . 
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E 




FIG. 5: The temperature dependence Tximp(T) for J = -0.2 and different t' = 0, 0.0005, 0.001, 0.002, 0.005, 0.01, 0.03, 0.05 
(lines from below to above) 
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FIG. 6: The temperature dependence Ximp (T) for different J and t' 



One can see from FiglS] that the minimum of x(^) occurs with decreasing A, its position being determined by 
the maximum of x^'^H^) Eq. ljlOp . This minimum is due to the strong energy dependence of the bare density of 
states, see Eas.([5|)- (ITT]) . Thus, besides the Kondo temperature Tfc, a second characteristic scale can occur in our 
problem, which is determined by A. Therefore we have to discuss the meaning of the Kondo temperature Tr- in more 
detail. This quantity comes from the expansion in J starting from high temperatures where impurity susceptibility 
x{T) obeys the Curie law. Tk is determined as a temperature scale where a considerable deviation from this law 
occurs. In the flat-band situation and for smooth p{E), x{T) decreases monotonously with increasing T, the universal 
behavior taking place. This picture is characterized by the ratio w = 4Tkx{^) ~ 0.41 0-I3l which relates high- and 
low-temperature scales Tk and x(0)- Such a behavior holds for Tk < A, but for small A values the universality is 
broken, so that the intermediate-temperature dependence x(2^) becomes complicated and w deviates from 0.41 (see 
Table III). 

The Wilson ratio R = (47r^/3)ximp(0)/7imp relating low-temperature susceptibility and linear specific heat 7imp is 
also presented in the Table H] It is important that even for rather small A the value of R is still close to 2. Therefore 
a Fermi-liquid behavior, characteristic for the Kondo problem is restored at low temperatures, except for the case 
of extremely small A where Ximp(O) and Cimp(O) can even become negative. 

Occurrence of negative values of Ximp and Snnp was demonstrated in Ref.[26| by a strict analytical consideration 
of the simple case J — —oo for the semielliptic density of states. This is a common property of systems with very 
narrow density of states peaks near the Fermi level. Of course, the quantities Xtot and Xband remain positive. 

In the standard situation of Kondo effect one has Tx{T) = 0{T), so that the impurity moment (note that Tx = 



(5f)t 



('S'z)band, see Appendix) is completely compensated by conduction electrons. At the same time, in the 
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TABLE I: Results of our NRG calculations for different t' and J: Tk (first line); Wilson ratio R (second line); the quantity 
w = 4rif x(0) (third hne; Wilson's value is w = 0.4128 ± 0.002) 



t' 


J = -0.1 


J = -0.15 


J= -0.2 


J= -0.3 




1.19-10"'' 


4.76-10"* 


0.00333 


0.0277 


0.05 


2.008 


2.008 


2.014 


1.99 




0.414 


0.416 


0.416 


0.415 




3.62-10"'' 


0.00103 


0.00613 


0.0450 


0.03 


1.995 


1.988 


2.005 


1.98 




0.413 


0.413 


0.414 


0.413 




3.45-10"" 


0.00436 


0.0184 


0.0930 


0.01 


1.998 


1.995 


1.991 


2.12 




0.415 


0.416 


0.425 


0.386 




7.49-10"* 


0.00726 


0.0288 


0.0989 


0.005 


1.999 


1.997 


2.069 


2.237 




0.416 


0.418 


0.447 


0.268736 




0.00112 


0.0102 


0.0320 


0.100007 


0.003 


2.002 


1.992 


1.998 


3.18 




0.417 


0.449 


0.392 


0.184 




0.00144 


0.0121 


0.0330 


0.1004 


0.002 


1.996 


1.97 


2.040 


-177" 




0.417 


0.44768 


0.337 


0.108 




0.00208 


0.0134 


0.0337 


0.1006 


0.001 


1.989 


2.03 


2.41 


0.343'' 




0.43 


0.38 


0.24 


-0.048 



" x(0) = 0.27, 7in,p = -0.02 
x(0) = -0.12, 7i„,p = -4.6 



pseudogap situation (low density of states near Ep, p{E) oc \E — £'i?|'",r > 0) the screening is incomplete, Tx{T) > 
[35| . Our problem describes the opposite situation: p{E — > 0) diverges and we have overcompensation: Tx{T) < 0. 

The singular case A — demonstrates essentially non-Fermi-liquid behavior (divergence of impurity magnetic 
susceptibility and specific heat at T — 0). This situation of singular density of states at the Fermi level is somewhat 
similar to the overscreened Kondo problem [36l [ where the number of scattering channel of conduction electrons 
n > 2S {S is the localized spin value); in this case \(T) demonstrates power-law behavior, but remains positive. The 
low-temperature behavior at A = can be fit as [2^ 



B 



'5imp = -TTZTTFTT^ 1.065,5 « 0.89. (27) 



(see Fig 12]), so that 



|in(r/i^)|*+i- (2^^ 



These asyniptotics are independent of J (cf. FigH]) since at sufficiently low temperatures any value of \J\ manyfold 
exceeds both the temperature and the width of the infinitely thin logarithmic peak in p{E). 

According to Ref.|37l]. the NFL behavior with C oc 1/ In^ (T/^/T) takes place in the case of underscreened {S > 1/2) 
Kondo problem [37| . 

Fig. [5] shows the temperature dependence of impurity specific heat Cimp = Ctot — Cband for different t' . As a 
rule, this dependence demonstrates two peaks. At not too small A, the high-temperature maximum occurs at the 
temperature, determined by the distance from VHS to Ep. This is owing to the non-monotonous dependence of 
X^°^(T), see Eqs.dTU]), ([131) • When decreasing temperature and passing this maximum, Cimp(-^) acquires a minimum 
and even can become negative. The low-temperature peak is owing to the Kondo effect and takes place in the standard 
flat- band situation too (see Ref. [111). One can see that its position corresponds roughly to the Kondo temperature. 
For small A < Tk, the order of positions of the maxima becomes interchanged. 

The corresponding magnetic entropy iSimp is shown in Fig. IHl One can see that this quantity tends to the value 
ln2 = ln(2S'-|- 1) at high temperatures and demonstrates the Kondo compensation at low temperatures. The behavior 
turns out to be non-monotonous due to the maximum in x^^H'^) (^^^ Eq. p3 
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FIG. 7: The fitting of thermodynamic properties at low temperatures for A = 



O 




FIG. 8: The temperature dependence of impurity specific heat for J = -0.2 and different t' = 0, 0.002, 0.005, 0.01, 0.03, 0.05 
(lines from below to above, if one sees the left-hand part of the figure); circles mark Tk, and triangles A 



CO 




FIG. 9: The temperature dependence of magnetic entropy for J = -0.2 and different t' = 0, 0.0005, 0.001, 0.002, 0.005, 0.01, 
0.03, 0.05 (lines from below to above) 
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IV. CONCLUSIONS 



In the general problem of magnetism of metals and alloys, the peaks in the bare density of states (which are usually 
owing to Van Hove singularities) near the Fermi level play a crucial role. The correlation effects can result in a 
considerable qualitative and quantitative modification of the corresponding temperature dependences [28|. They can 
be important, e.g., for transition-metal alloys. 

The considered t — t' Kondo problem which includes VHS is a non-trivial example of influence of density-of-states 
peaks on electron properties being combined with correlation effects [28l Isij ] . Our treatment gives an example of exact 
numerical solution of such a problem. The resulting temperature dependences of thermodynamic properties include 
both one-particle effects, connected with VHS, and many-electron Kondo features. Therefore the overall picture is 
rather complicated. At low temperatures, the Fermi-liquid behavior is restored, except for the case of very small t' . 
In the latter case, a non-Fermi-liquid behavior takes place which should be studied by more advanced methods. 

It would be also of interest to perform similar calculations for the Kondo lattice problem, e.g., in some "mean- field" 
approximation. A "poor man scaling" approach was applied to this problem in Ref . (ssj . 

The research described was supported in part by the Program "Quantum Physics of Condensed Matter" from 
Presidium of Russian Academy of Sciences. The author are grateful to A.O. Anokhin and A. A. Katanin for useful 
discussions. 



Appendix. Numerical renormalization group approach for the singular density of states 



Here we discuss some important details of the numerical renormalization group (NRG) method @, ^ as applied 
to our problem of the singular density of states. 



Construction of the Wilson chain 



Following to Wilson Q| we use a unitary transformation to pass from the operators Ck to the operators /«. Then 
the impurity model with a Hamiltonian of type ([ij is reduced to a semiinfinite chain (Fig llOp with a Hamiltonian of 
type: 



H 



sd 



OO 



(29) 



The renormalization group procedure starts from the solution of the isolated-impurity problem (sites "imp" and eo 
in Fig. nop . At the first step, we add a first conducting electronic site ei, and construct and diagonalize a Hamiltonian 



J 70 71 72 
• '-O^ '-O^ i-O 

imp eo ci <^2 £3 



FIG. 10: Representation of the Kondo model in the form of a semiinfinite Wilson chain 



matrix on this Hilbert space (with a fourfold higher dimensionality). This procedure is multiply repeated. However, 
since the dimensionality of Hilbert space grows as 4^ {N is the order number of an iteration), it is impossible to 
store all the eigenstates during the calculation. Therefore, it is necessary to retain after each iteration only the states 
with the lowest energies. If we restrict ourselves to a certain maximum number of stored states (determined by the 
computational possibilities), it is necessary, starting from a certain iteration, to leave of the order of 1/4 of states at 
each step. 

Unfortunately, direct application of this scheme fails, since the disturbance introduced by the elimination of the 
high-lying states proves to be too large. Wilson found a method of overcoming this difficulty. This reduces to the 
artificial introduction of an exponential suppression of matrix elements 7„ , which decreases the coupling between the 
retained and eliminated states, thereby decreasing the infiuence of the eliminated states. To this end, Wilson [1] used a 
logarithmic discretization of the conduction band, i.e., replacement in ^ of an energy range G [ryDA"™, riDA~™-~^^], 



11 



rj = 1, —1, m = 1, 2, 3, ... by a single level with an energy equal to the average energy of this interval (D is the 
half- width of the conduction band, A > 1). This results in a change in the density of states: 

) , (30) 

r/,m 

where 

As a result, the jumps will have the required decay, 7„ cx A^"/^; for a flat band, Wilson obtained analytically 

i:>(i + A-i) (1- A-"-i) 



7« 



2Vl - A-2"- Vl - A-2«- 



A-"/2,e„ = 0. (31) 



In more complicated situations, the construction of the Wilson chain must be performed numerically. Usually (see, 
e.g., Ref.[4l|) this is performed in spirit of the initial work [2], by numerical reproducing Wilson's analytical flat-band 
procedure for an arbitrary density of states. However, there exists another way which yields equivalent results, but 
seems to be more natural. To bring the Hamiltonian matrix 

to the tridiagonal form, one can use the Lanczos tridiagonalization algorithm [4^ which is just adapted namely 
for this problem. For the model (H]) this method was described in Ref.[J| where a Wilson chain for non-discretized 
semielliptic density of states was also analytically constructed. We perform this procedure numerically for an arbitrary 
logarithmically discretized density of states. Starting from the vector |0) = /oo.|vac) (where /q^ = J^n m ^vm.cljn,a) 
we generate a new basis |0), |1), |2)... for the conduction electron states by Schmidt orthogonalization: 

|1) -(i/.|0)-|0)(0|i7.|0)) , 

70 

|2) = -iH,\l) - \1){1\H„\1) |0)(0|iJ.|l)) , (33) 

71 

\n+l) = —{H,\n)-\7i,){n\H„\7i,)~\n-l){72-l\H,\n)), 

In 

where each 7„ is chosen to normalize \n + \). One can see that {n'\Hrj\n) = for n' = 0, 1, 2...n — 2. This means 
that Her is tridiagonal in the new basis. The off-diagonal elements are {n + l\Hcy\n) = 7„. Defining e„ = {n\Hcr\n) we 
consecutively obtain the coefficients 7„ and e„ for ((29| . 



Calculation of thermodynamical averages 

In our calculations we put A — 1.5 and store 10^ states per iteration. 

Diagonalizing the Hamiltonian (|29|) for a given chain length N yields a set of eigenvalues. As indicated in Ref.[2l, 
because of retaining only part of the energy spectrum at the A^-th step of the NRG procedure, thermodynamic averages 
should be calculated at a temperature that depends on A: T/v = A^^/^Tq, where the starting temperature Tq is chosen 
more or less arbitrarily. However, it should be neither too large (then the contributions of abandoned high-energy 
states becomes important), nor too small (then the discreteness of the energy spectrum becomes appreciable). 

The total entropy and specific heat read [i^ 

Stot = {H)tot/T + \nZtot , 

Ctot - [(i/')tot - {H)l,] IT^ , (34) 

where Z is partition function. Since the total spin commutes with the Hamiltonian H , each eigenvalue is characterized 
by a well-defined spin projection Sz^tot, therefore the quantities like {SDtot can be calculated straightforwardly. On 
differentiating {Sz)tot with respect to magnetic field one obtains [2] 



Txtot(T) = [{S!},u - {Sz)U 



(35) 
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T'Xband(T'), iSband and Cband are Calculated in a similar way, and the corresponding impurity contributions are obtained 
by subtracting them from ([M)) - ([55t . 

Because of finite length of the chain, the thermodynamic quantities like TVximp(7A') and Cimp(7V) demonstrate 
even-odd oscillations depending on T^r which have nearly constant amplitude. Therefore the amplitude of the oscil- 
lations in x(T) and C{T)/T = ^{T) increases strongly with lowering T. 

To suppress the oscillations in x(T), we used smoothing according to Euler If there is a certain oscillating 

sequence we introduce a new sequence An^ whose members are equal to averages of the adjacent members of 
the initial sequence: A^^ = [An -\- A„+i)/2. If necessary, the procedure is repeated: An'' = [An'' + A^^'l^)/2. In 
particular, this was made in the calculation of Ximp- By designating xn = Ximp(Ijv), we obtain 

X^^\\/TnTn+i) = ^XN + \XN+1 
X^'^ {Tn) = \XN-1 + \XN + \XN+1 (36) 

X^^\\/TnTn+i) = \XN-1 + IXN + |XW+1 + \XN+2 

A similar problem which occurs at calculating the slope of specific heat 7imp is solved by the same way. The method was 
tested for the flat-band case to obtain the values R = 2.008, 2.016 and w = 0.416, 0.417 for J = —0.1, —0.2 respectively 
(cf. Table 1). Our method differs from [5] by that we calculate x and 7 directly rather than by constructing an effective 
Hamiltonian explicitly near the fixed point J ~ —00. Although resulting in a slight decrease of accuracy, such an 
approach can be applied more widely, in particular to obtain a NFL behavior (see Figs|4H6]). 
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